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The classical transitions between topologically distinct vacua in a Sf7(2)— Higgs 
model, using a Higgs field of mass approximately 120 GeV, is examined to probe 
the crossover region between the symmetric and broken phase. For the volumes 
used, this crossover is approximately 10 GeV wide. 



1 Introduction 

lis clear now from a number of numerical simulations in 3 and 4 dimensions 
that the first order phase transition from the broken to symmetric phase 
of the electroweak sector of the Standard Model ends for a Higgs mass of 
approximately 80 GeV. For larger Higgs masses, there is analytic crossover 
or very possibly a phase transition of third or lower order (whether numerical 
studies could ever distinguish between these two scenarios remains to be seen) . 
While the present quoted lower bound of 65 GeV for the Higgs jpass U and 
an estimate for the Higgs mass from electroweak precision dataci is rriH ~ 
145^77'' GeV is still consistent with an explanation of baryogenesis via this 
phase transition, the amount of parameter space is rapidly shrinking. 

It seems logical to ask what is occuring to the real time process of baryon 
number violation for larger Higgs masses. Does the crossover become increas- 
ingly broad as the Higgs mass is increased or does it remain reasonably sharp 
(changing from the symmetric to broken phase over a few GeV) ? An inter- 
esting scenario, discussed elsewhere in these proceedings El, suggests that for 
a different cosmology at electroweak temperatures one can still generate sufh- 
cient levels of baryon number violation, even with crossover. At the very least, 
understanding baryon number violation for this model will be important if we 
wish to employ a similar mechanism for more complicated field theories such 
as SUSY or Grand Unified Theories. 



2 Real time transitions 

Many of the details of this calculation have been ijavered extensively in these 
proceedings and so I will not go into detail hereB'ElEl. Suffice it to say that we 
wish to examine the real time process of the rate of change of the Chern-Simons 
number, which is related to baryon number violation. While a full quantum 
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field theory calculation is almost impossible numerically, a classical treatment 
of the problem may give a reasonable estimate for the rate. 



2.1 Classical analysis 

The following lattice SU (2)-Higgs Hamiltonian was constructed 



^fl-lrrf/p] -i^(A,$)t(A,cf>) (1) 
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The sums are defined over spatial sites only and the sum corresponds 
to a sum over spatial plaquettes of the gauge field Ui. The first two terms in 
the Hamiltonian are discretized versions of F^jF^j and {Di^y Di^ respectively. 
The fields Ei and tt are the conjugate momenta to the coordinate fields Ui and 
$. The lattice spacing, a has been taken to be 1. Formally, Eqn.(|^) results from 
the dimensional reduction of the 4-dimensional S'[/(2)-Higgs Lagrangian with 
kinetic terms reintroduced. One would therefore expect in the hamiltonian 
a Debye mass term and a renormalisation of the kinetic terms. As noted by 
Moore and Turoklj, the Debye mass term is effectively introduced by the U.V. 
cutoff of the lattice. The renormalisation of the kinetic terms is expressed in 
the coefficient ze and z^r. It is believed that they take the form 1 + 0{g^) and 
for this calculation they are assumed to be 1 . The time step factor At is taken 
to be 0.05, which has been used in similar simulations. 

A set of initial configurations are generated which satisfy Gauss constraints 
and then evolved classically using a leapfrog algorithm. Instead of evaluating 
the Chern-Simons number directly, which have a number of inherent problems, 
it is indirectly measured using a slave-field method, proposed by Moore and 
Turok Ej. This method has several advantages, being computationally quite 
cheap and does not require a finite renormalisation. 

One can set the Higgs mass by fixing the self-coupling term Xe since at 
tree level 

« 2Ai , (2) 
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for niH ~ 120 GeV, Xe was set to 1.125. In order to get an estimate of the 
temperatures, the relationships between the 4-dimensional and 3-dimensional 
parameters can be usedO to determine the ratio mn/T as a function of the 
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Figure 1: DifTusion rate of winding number at /3 = 7.1 and /3 = 7.6. Note the difference in 
the vertical scales. 

remaining bare parameters, /3 and Mm- For mn ~ 120 GeV the temperature 
of crossover, Tcross is around 200 GeVB. Hence, Mho was set to -0.596 so that 
aTcross corresponded to /3 w 7.5 (which is inversely proportional to aT). 



2.2 Results 

For a lattice size of 24^, (3 was varied from 6.8 to 7.8. The number of initial 
configurations and number of sweeps are shown in table |l|. As can be seen from 
figure 1^, one could determine the diffusion coefficient directly from the data by 
evaluating the slope. However, in order to maximise the sample size, a cosine 
transform was applied Ii3 to the data, and the first 200 coefficients evaluated. 
The diffusive contribution to the coefficient will be proportional to l/m?, where 
m is the number of the coefficient. An example of the resulting coefficients is 
shown in figure |^. The error for each coefficient is evaluated via jack-knife. 
Correlations between different coefficients appear to drop significantly as the 
statistics is improved with the percentage of correlation coefficients greater 
than 0.9 dropping to less than 5% for the larger data sets. The coefficients 
were averaged in bins of 50, which reduces the error significantly. The errors 
were summed quadratically. (One could also evaluate the error for the binned 
data by treating the central values of the coefficients as independent points 
and estimating the error as a standard deviation. This approach varied the 
statistical error by a factor of 0.5 to 2. Ideally, one should perform a correlated 
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p 


^configurations 


#tinie steps 


K 


6.8 


16 


7000 


0.973 ± 0.057 ± 0.161 


6.9 


16 


7000 


0.984 ± 0.061 ± 0.107 


7.0 


16 


10000 


1.013 ± 0.083 ± 0.187 


7.1 


8 


20000 


0.770 ± 0.059 


7.2 


8 


10000 


0.388 ± 0.034 ± 0.038 


7.3 


8 


10000 


0.092 ± 0.008 


7.4 


8 


6000 


0.018 ± 0.002 ± 0.003 


7.5 


8 


9000 


0.00787 ± 0.00068 


7.6 


16 


6000 


0.00088 ± 0.00006 


7.7 


8 


6000 


0.00130 ± 0.00016 


7.8 


8 


9000 


0.00055 ± 0.00005 



Table 1: Numbers of configurations and the time steps iterated forward for each /3 . The 
second error quoted for k is the difference between the first and second bins of 50 coefficients 
when the difference was greater than 1 standard deviation. 

fit through these points, however the sample size is not large enough to do this.) 
For some values of /3, differences between different regions of binning varied by 
more than one standard deviation. A conservative estimate of the systematic 
error due to fitting took this into account. 



3 Conclusions 

From figure || one can see that the diffusion rate vanishes for (3 > 7.4. The 
diffusion rate is decreasing from jS = 7.0. This corresponds to a temperature 
variation of around 5%, which, assuming the crossover temperature is around 
200 GeV implies a variation of about 10 GeV. The ra±e Jql, (3 = 7.0 is of the 
same order as similar simulations done for niH ~ mw E3l3'c2l. The calculation 
of K for P < 7.0 is consistent with k{(3 = 7.0); however the Moore- Turok def- 
inition of the winding number will very often change by more than 1 between 
measurements which introduces a systematic error. The next step is to adjust 
Mho so that aTcross is much smaller and thereby larger /3's (smaller lattice 
spacings) are required. A study of the effect of the finite volume is also nec- 
essary. Ultimately this should be repeated for a larger Higgs mass (150 GeV, 
for example) to determine if the width is increasing with the Higgs mass. 
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Figure 2: The mean square of the cosine coefficients at /3 = 7.4 and the average of the 
coefficients with a bin size of 50. 



— * — $ — $ — $ — ' — I 
7.5 6.0 



Figure 3: k as a function of /3 using a lattice. 



5 



Acknowledgments 



I am indebted to Guy Moore and Neil Turok for making their code available 
for this simulation. I thank Anne Davis, Ron Horgan, Alex Krasnitz and Jan 
Smit for useful conversations. This research is supported by the Leverhulme 
Foundation. The numerical work was carried out on a Silicon Graphics/Cray 
Origin 2000 run by the UK Computational Cosmology Consortium and an 
Hitachi SR2201 at the High Performance Computing Facility at the University 
of Cambridge. Support for these computational facilities has been provided 
by Silicon Graphics/ Cray, Hitachi, HEFCE, PPARC and the University of 
Cambridge. 



References 

1. F. Csikor et al, Nucl. Phys. (Proc Suppl.) 53, 612 (1997). 

2. K. Kajantie et al, Phys. Rev. Lett. 77, 2287 (1997). 

3. Particle Data Group, R.M. Barnett et al, Phys. Rev. D 54, 1 (1996). 

4. J. Ellis, G.L. Fogh and E. Lisi, Phys. Lett. B 389, 321 (1996). 

5. M. Joyce, these proceedings. 

6. G. Moore, these proceedings. 

7. A. Krasnitz, these proceedings. 

8. J. Smit, these proceedings. 

9. G. Moore and N. Turok, Phys. Rev. D 55, 6538 (1997)^ 



10. G. Moore and N. Turok, PUPT-1681, DAMTP-97-17, [licp-ph/9703266 
(1997) 

11. K. Kajantie et al, Nucl. Phys. B 458, 90 (1996). 

12. W. Tang and J. Smit, Nucl. Phys. B 482, 265 (1996). 



6 



